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Abstract 

Following the ideas of L. Carlitz we introduce a generalization of the Bernoulli and Eulerian 
polynomials of higher order to vectorial index and argument. These polynomials are used for 
computation of the vector partition function W(s,T)), i.e., a number of integer solutions to a 
linear system x > 0, Dx = s. It is shown that W(s, D) can be expressed through the vector 
Bernoulli polynomials of higher order. 

1 Introduction 

The history of the Bernoulli polynomials B^{x) counts more than 250 years, as L. Euler first studied 
them for arbitrary values of the argument. He introduced the Euler polynomials Ek(x), and relation 
between these functions was established in the end of nineteenth century by L. Saalschutz. The 
generalization of the Euler polynomials to the so-called Eulerian polynomials Hk(x,p) was made by 
Frobenius [7] several years later. 

N. Norlund [9] introduced the Bernoulli £?[ (a;|d m ) and Euler polynomials E^ n \x\d m ) of higher 
order adding m parameters. Similar extension for the Eulerian polynomials H^ n \x, p m \d m ) was 
made by L. Carlitz in [3]. 

The Bernoulli and Eulerian polynomials of higher order appear to be useful for the description 
of the restricted partition function W(s,d m ), which is a number of integer nonnegative solutions 
of Diophantine equation d m • x = s (see The author showed in [10] that W(s,d m ) may be 

written as a finite sum of the Bernoulli polynomial of higher order multiplied by prime circulator 
introduced by A. Cayley (see [6]). The short derivation of this result is given in Section [3j 

Carlitz suggested in [5 J another extension of the Bernoulli -Bk( x ) and Eulerian polynomials 
i?k( x ) p) to vectorial arguments and indices. 

In this work we combine both abovementioned approaches to introduce a new class of polynomi- 
als, which we call vector Bernoulli S^' m ^(x|D) and Eulerian H^' m ' (pc, p\T)) polynomials of higher 
order. We define the generating functions for these polynomials, find the recursive and symmetry 
relations, and find a new relation between these polynomials. Using the vector polynomials of higher 
order in Section |4] we find an explicit formula for the vector partition function W(s, D), defined as 
a number of integer nonnegative solutions to the linear system x > 0, Dx = s, where D denotes 
a nonnegative non-degenerate integer matrix. It appears that similar to the scalar case the vector 
partition function may be written as a finite sum of the vector Bernoulli polynomials of higher order 
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multiplied by prime circulators of vector index and argument. The solution gives the vector parti- 
tion function in every chamber of the system, as well as it determines the shape of each individual 
chamber. 



2 Bernoulli and Eulerian Polynomials and Their Generalizations 

In the Section 12.11 the definition and main properties of the (regular) Bernoulli B k {x) and Eule- 
rian H k (x,p) polynomials are presented. Their generalization to the polynomials of higher orders 
B { k m) (x\d m ) and H { ™\x lP m \d m ), recursive and symmetry properties and some useful formulas are 
presented in Section 12.21 Extension of the regular polynomials to vectorial arguments and indices 
-Bk(x) and H^x, p) described in the Section 12.31 Finally, in the Section 12.41 combining both ex- 
tensions we introduce a new class of the polynomials f?^ (x|D) and (x, p|D) and consider 
their properties. 



2.1 Bernoulli and Eulerian Polynomials 

Define the Bernoulli polynomial B k {x) through the generating function 



k=0 



The Bernoulli numbers are defined as B k = B k (0). It is known that the multiplication formula holds 
for these polynomials: 

p-l 



Vi? fc (x + -) =p 1 - k B k (px); 
J=o V PJ 



this relation can be used as an alternative definition [8] of the Bernoulli polynomials. Another 
remarkable property of the Bernoulli polynomials is that they can be presented as a symbolic power 



B k (x) = (B + x) k = ^( k )B 



x k ~ 



where after the expansion the powers of B are transformed into the indices. 

The Eulerian polynomials were considered by Frobenius in [7], their generating function reads 

H k=0 

For p = — 1 the Eulerian polynomials reduce to the Euler polynomials E k {x). The former also 
admits the symbolic form 

H k (x,p) = (H(p)+x) k , 

where H k (p) = H k (0,p) are the Eulerian numbers. For pP = 1 we have an additional relation [3] 
between the Bernoulli and Eulerian polynomials 

p k-iJ2p- r B k (x+-) =J^H k _ 1 (px,p). (1) 
f^o V PJ 1 ~P 
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2.2 Bernoulli and Eulerian Polynomials of Higher Order 

Norlund [9] introduced the Bernoulli polynomials of higher order l?( m ) (x|d m ) adding m parameters 
d m = {di, d>2, ■ ■ ■ , d m } by means of the generating function: 



m j 00 + k 



i=l fc=0 

The regular Bernoulli polynomials are expressed through the relation 



B<£\x\d) = d k B k (x). 
Norlund also found the multiplication formulas 

VM m) (x + -c7(d m )|d m =p m - fc ^ m) (px|d m ), a(d m ) = Vd 4 , 

Pi— I in 

J2 Bt\* + r • & m \{vid U ■ ■ ■ ,Prnd m }) = 7r(p m )B^\x\d m ), vr(p m ) = \{p u 

r;=0 i=l 

where r • d m = Y^iLi denotes the scalar product. Norlund gave the symbolic expression 

m 

B { ™\x\d m ) = (x + J2di i B) k , 



i=l 



where again after the expansion powers of l B are converted into the indices of the Bernoulli numbers. 
This relation implies a binomial formula: 



(2) 



k=0 

The following symmetry relation is valid 

B { k m \x\d m ) = {-l) k B { k m \-x - <r(d m )|d m ), (3) 

and the recursive relation holds for the Bernoulli polynomials 

B { k m \x + d m \d m ) - B ( k m \x\d m ) = kdmB^ix]^- 1 ). (4) 

The Eulerian polynomials of higher order H^ m \x, p m \d m ) were considered by Carlitz in [3]. 
Their generating function reads 

00 



i=l Pi fc=0 



and they admit the following symbolic notation: 

m 

H^\x,p m \d m ) = (x + J2diH( Pi )) k . 



i=l 



When py = 1, pj 7^ 1, j = 1, 2, . . . , m we have the analog of ([!]) 

— — £ PpB<ri(x + r • d TO |{pidl,P2d2, • • • ,Pmd m }) 
\" ' »■ — n 



fc! / TT ft^i \ rr(m) 



ci(^i dm )- (5) 



(k — m)\ I , 

v \j=i 

2.3 Vector Bernoulli and Eulerian Polynomials 

Carlitz suggested another extension of the Bernoulli and Eulerian polynomials. He introduced 
[5] the Bernoulli polynomials -Bk( x ) of the vector argument x = {x\, . . . , xi) and vector index 
k = {k\, . . . , ki} through the generating function 

e Z^=i H — \ k K. 



where the summation is over all k\, ■ ■ . , k\ and x • t = Yli=i x i^i is a scalar product. Here and below 
we use the notation 



i=i i=i v 7 i=i y 1 



The multiplication theorem is established in the form 

p-l 



V B k ( x + -) = pi-MflkCpx), |k| = V A*. 



The Eulerian polynomials H^x, p) of the vector argument and index are defined as 



xt 



1 00 t' 4 

e2-<=i*._p k=Q 



and the analog of ([I]) for pP = 1 reads: 

p|k|-l £ p -r fik ^ + ^_ J- kiH^v (px, P ), 

where 1® is vector with the components 1„ = <5j m , and 5j m is the Kroneker delta. 

Here and below we will use the notion of the vector Bernoulli B^(x) and Eulerian H^(x, p) 
polynomials remembering that they are scalar functions of their arguments. 

2.4 Vector Bernoulli and Eulerian Polynomials of higher order 

In this work we make one more step in generalization of the Bernoulli and Eulerian polynomials which 
naturally arises from the approaches discussed above. We consider vector Bernoulli B^ ,m \x\'D m ) 
and Eulerian (x, p m |D m ) polynomials of higher order, where D m denotes (I x m) matrix. 



4 



The vector Bernoulli polynomials of higher order are defined by means of the generating function: 

m v z t n - - + k 

2_/i=i L i 1Jl 



/i^^-l'V (X|D } k! 



Denoting the columns of the matrix D m as {ci, . . . , c m }, we can rewrite the above definition as 

m , k 



3=1 k 

The analog of the recursive relation reads 



t n^ T = E^ m) (x|D-)^ [ . (6) 



i=l 

One can check the validity of an analog of symmetry relation Q 

m 

^ ,m) (x|D m ) = (-l)l k l^' m) (-x-cr(D m )|D m ), cr(D m ) = J^c 

Using we obtain the vector generalization of the binomial formula ([2]) 

Q^' mi) (x|D mi )B^ 2) (y|D m2 ) = ^' mi+m2) (x + y|D mi 



3=1 



m,2 N 



(7) 



k 

The vector Eulerian polynomials of higher order have the generating function of the form: 

m i j-k 

3=1 Pj k 

where p m = {p\, . . . ,p m }- Consider a relation between the vector Bernoulli and Eulerian polyno- 



mials of higher order similar to ©. For pj such that p P : j = 1 we have the following expression: 



Pj-l m k 

Yl E pJ rjB k( x + 5Zrj c il{pi c i> • • • ,p m c m })- 



k! 

i=i 



./ i '•, 11 i=i 



in 



(m \ , k 

n^(c,-t)j£/4nx,pHD m 4 



*(P m ) T^r(^-t) > ^r ; (x,p m |D m )-. (8) 



3 Restricted Partition Function 



The restricted partition function W(s, d m ) = W(s, {d±, cfo, ■ ■ ■ , d m }) is a number of partitions of an 
integer s into positive integers {d±, cfo, . . . , d m }, each not greater than s. The generating function 
for W(s,d m ) has a form 



III LXU 



i=l s=0 

where W(s, d m ) satisfies the basic recursive relation 

W(s, d m ) - W(s - d m , d m ) = W(s, d™" 1 ) . (9) 

Note a similarity of Q to the recursive relation ([¥]) for the Bernoulli polynomials of higher order. 
Sylvester found a symmetry property of the partition function: 

m 

W(s,d m ) = {-l) m - 1 W(-s-a(d m ),d m ), a(d m ) = Y,di- (10) 

i=l 

He proved [12] a statement about splitting of the partition function into periodic and non-periodic 
parts and showed that the restricted partition function may be presented as a sum of " waves" , which 
we call the Sylvester waves 

W( S ,d m ) = ^2w j (s,d m ), (11) 

i=i 

where summation runs over all distinct factors of the elements of the set d m . The wave Wj(s, d m ) is 
a quasipolynomial in s closely related to prime roots p~ of unity. The wave Wj(s, d m ) is a coefficient 
of t~ l in the series expansion in ascending powers of t of the generator 

F™M) = E 7^-, (12) 

The summation is made over all prime roots of unity pj = exp(2irin/ j) for n relatively prime to 
j (including unity) and smaller than j. It is easy to check by straightforward calculation that the 
recursive relation 

F™(s,t) - F™(s - d m ,t) = F™-\s,t) (13) 

holds for any generator FJ n (s,t), impliyng the validity of (J9j) for each Sylverster wave Wj(s,d m ). 
The generator satisfies the following symmetry property 

*f (s, t) = (-ir F p(-s - a(d m ), ~t), 

which implies the validity of (PT0|) for the residue of ^^(s,*). 

In [llj the explicit expression for the Sylvester wave of the arbitrary period is given through 
the Bernoulli and Eulerian polynomials of higher order. Using ([5]) it was shown in [10] that it is 
possible to express the Sylvester wave as a finite sum of the Bernoulli polynomials of higher order 
only. Here we present a short derivation of this result. 

Assuming that the vector d m has luj components divisible by j, sort the elements of d m in such 
way that the j-divisible integers come first. The generator (|12p can be written as a product 

Fj(s,t) 



~ n£i(i - x n^. + i(i - P fe- dit ) 
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Consider a modified generator F^^s^t) = TT(d m )t m FJ n (s,t), for which using the notation p d " 
{p dl ,p d2 ,..., p dm } we have 



Pj i=i ■ 1=^+1 e * 



t i m -i _ ai j a; 

-S-<r(d"' "J)TT gif Js+a(d m ))t T T O 



X 



t m-w J V^ -«-«r(d J )TT ^ e ( S +^(d m ))t TT 

p j 11 e d 4 t _ i 11 di t _ n <n i _ 

Pj i=l i=wj+l e 1 

t m-^ V- - S -<x(d— i) . y B to) (0 | d «*)* (14) 
^— ' -' ni! 

Pj ni=0 

II r\ • E + '(o,/T^ i d "-" j )b- 

i=u>j+l 1 n 2 =0 2 ' 



Now we employ ([5]), noting that for lo + 1 < i < m all pi = j , to obtain 

m j di oo n 

n ^■y. h ^ 3 \ s +<^pT^\^)- 

- LJ - 1 — n. — no\ 

i=Uj+l 1 n 2 =0 A 

- E (m _^ +n2)!jm -. J E ,7*'%ZU> + «*n + ■ ■ w 

n 2 =0 v r i= 

Inserting the above expression into (|14p and using ([2]) we obtain 

3-1 



^ m ( S ,t) = J - (m -^ ) EE'°7" r ' d " J ^ (dm) E 5 n m) ( s+cj ( dm )+ r - dm "" 3 i d i 1 )^' ( 15 ) 

Pj r i= n 



where we use a shorthand notation d™ for j -modified set of summands defined as union of subset 
d w J of summands divisible by j and the remaining part d m ~ u i multiplied by j 

df = d u * Ujd m ^ = {d h . . .,d Wj JtLj+u . . . ,jdm}, 

so that d™ is divisible by j. From (fT5|) the expression for Fj n (s,t) follows: 

f ;(M) %(d i m -, EE^r j -^^4™) (s+C x(d™) + r.d^id r )i- r . (16) 

V ^ pj r i= n 

Setting in n = m — 1 we arrive at the expression for the Sylvester wave: 

1 j ~ l 

WJs,d m ) = - L - Y B ( ™\{s + o-{d m ) + r ■ d m ~^\d™)Y p~: s - r - Am ^^ Am \ 

•? v ' ' (m - 1)! 7r(d m ) j m_a; J m ~ l 3 ; Z^O 

ri=0 Pj 

Introduce a notation 

pj 

for the prime radical circulator (see [6]). For prime j it is given by 



<Ki), s = (modj), 
(modj), 
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where cj>(j) and fj,(j) denote Euler totient and Mobius functions. Considering j as a product of 
powers of distinct prime factors 

k 

one may easily check that for integer values of s 

k \Pk / 

where ^fc(s) = for non-integer values of s. 
Noting that for j-modified set d™ we have 

m 

7r(d™)=j m -^vr(d m ), a(df) = a(d m ) + (j - 1) £ 
and using the prime circulator notation we can write the Sylvester wave in a form 

wj{;* m ) = (m _ ^ ff(d?>) E B ^( s + CT ( d D - r • d^|d-)^( S - r • d— o- 

The polynomial part of the partition function corresponds to j = 1 and equals to 

W 1 (s,d-) = (m _ 1) 1 , 7r(dm) ^ m - ) i(^ + ^(d m )|d-). (17) 
The polynomial part for the j-modified set d™ for j > 1 reads 

^(-.dy) = (m _ 1) ! 7r(dr) - B i m - ) i( g + ^df)|d?)- (18) 

Thus, the Sylvester wave for j > 1 can be written as a linear superposition of the polynomial part 
of the j-modified set multiplied by the corresponding prime circulator: 

j'-i 

Wj(s, d m ) = ^ Wi(a - r ■ d m "^, d™)^(s - r • d m ~^)- (19) 
n=o 

It is easy to see that for j = 1 one has d™ = d m , ^i(s) = 1 and u>\ = m, so that summation signs 
disappear, and (fT9|) reduces to ([T7]l . 



4 Restricted Vector Partition Function 

Consider a function W(s, D m ) counting the number of integer nonnegative solutions x > to the 
linear system D m x = s, where D m is a nonnegative integer (I x m) matrix. The function W(s, D m ) 
is called vector partition function as it is natural generalization of the restricted partition function 
to the vector argument. 

The generating function for the vector partition function reads 

m -. 

n 7— ^ = e w & Dm ) ts = e w & {*>■■■> c ™» ts . ( 2 °) 

i=l s s 
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where c.j denotes the i-th. column of the matrix D m . It is easy to see that the vector analog of the 
recursive relation ([9]) holds 

W(s, D m ) - W(s - c m , D m ) = W(s, D'"- 1 ). (21) 
The symmetry property established in pQ is just a vector generalization of (jlOp : 

m 

W(s,-D m ) = (-l) m - rankDm Ty(-s-cr(D m ),D m ), cr(D m )=^c i . (22) 

i=i 

Similarly to the Sylvester splitting theorem (jlip we write the vector partition function as a sum 

of vector Sylvester waves 

W(s,D m )=^Wj(s,D m )= £ ^j(s,D m ), (23) 

where j denotes /-dimensional vector (ji , . . . , ji ) . The summation for each runs over all distinct 
factors of the elements of the matrix D m . The vector analog of the generator (|12p is written in the 
form (see [3]) 



where 



F- m fs t) — \ — 

pr=n^ H - ( 24 ) 



fc=i 

It satisfies the relation similar to (|13p 

J Fj m (s,t)- J Fj m (s-c m ,t)=F j m - 1 (s,t) 
for any generator i^ m (s,t). The following relation aslo holds 

Fj m (s,t) = (-l) m i^ m (-s-cr(D m ),-t). 

The multidimensional residue of the generator gives Wj(s, D m ) as the coefficient of t" 1 = Y[ k= i l k 
and each vector wave satisfies the relations ()2ip and (|22p . 

Let an equality pj l = 1 holds for oj- } columns of the matrix D m , and sort the matrix in such 
way that these Uj columns come first, so that 1 < i < Wj. Introduce a homogeneous polynomial of 
degree m 

m 

P m (t,D m ) = [](c,,-t)= £ C N (D m )t N , 

i=l |N|=m 

and construct the modified generator 

i^(s,t) = P m (t,D m )i^(s,t), 

which can be written as 

j 777, -i Cj 

^(s, t) = Eprn^x^^.n c, " t 



p. 1 i - p. 1 

p Jn i=l i=u>j+l J J 



J ^ — ' ri i ! 

in m 

II ttt- E ^r^ j) (s + -(D m ), pf-^ id— J -)i T , 
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where elements of the vector p m W J are pj* 7^ 1 given by (|24p for uj- } + 1 < i < m. Using the relation 
([8]) with pi such that p^ lCl = 1 we can write 

ft E ^^(s + er(D»), p? |D— i)£ 

1 (; - \ m t" 2 

= -rfom-aM E E PP^^ ^ ( S + E »'iCi|{Pa, J +lC Wj+ l,...,p m C m })— T . 
VH 7 n 2 r i= i=u>,H-l 



Employing the binomial formula ([7]) we obtain i^ m (s,t) in the form 

Pp( S ,t) ~EE B " m) ( s + ff ( Dra )+ r ' Dra "i D r)/Er ,D ' ; 

^ p ; n n=0 " p J4 

where Df 1, is the i-modified matrix of the form 

D]" = {Ci, . . . ,^,^.+10^+1, . . . ,p m Cm}. 

Noting that for the j-modified matrix D™ 

m 

P m (t,Df) = ^(p m -^)P m (t,D m ), ( r(D J m )=<T(D m )+ J2 fo- 1 )^ 

i=u)j+l 

and using the vector prime circulator notation 

/ 

fc=i 



we write for the generator i^ m (s,t) 

P; m (s,t) = p-^Df) E S n ( s + °"( D r) - r • D m -^|D] n )^ j (s - r • D m -^) — . 



n! 

n r;=0 



In the above expression only terms with |n| = m — I contribute to the vector Sylvester wave 
Wj(s, D m ), which is found also as a constant term of -Fj m (s,t)t 1 , equal to a fraction with both 
numerator and denominator being homogeneous polynomials of degree m 



W j (s,D m ) = limP m 1 (t,Df) £ ^P^VlDf^s')-^, s' = s + <x(Df)-r.D' 



It is convenient to write each vector Sylvester wave Wj(s, D m ) as a sum of quasipolynomials, 
which we call partial vector Sylvester waves: 

VKj(s,D m )= Wj n (s,D m ). (25) 

\n\=m—l 

The number of different partial waves is found as a number of representations of an integer m — I 
as a sum of I nonnegative integers, i.e., it is equal to 



W(m-l,l l 



m — 1 
m — I 
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The partial wave W^ n (s,D m ) is found as (|n| = m — l) 

^(s,D m ) = limP-^t.DH V ^' m) (s + cr(Df) - r-D m -^|D^)^i(s-r -D m ^)— -. 

J J / ^ J J « j--^J 



wi 



Using a parametrization t = at we compute the above limit as 

Pi-l n+l 

l^ n (s,D m ) = p-^Df) ^ J B^' m) (s + o-(D] n )-r-D m - w -i|Df)^j(s-r-D m ^) r 

n=o n ' 

Pi— l 

= C n (a,Df) ^ P^' m) (s + cr(D] n ) - r • D m ^ |Df )*j(s - r • D m ^), 

n=o 

where 

a n+1 

CJa, Df 1 ) = -; 7 r, n = m - Z. 

We define the polynomial part of the vector partition function for j = 1, the corresponding 
partial polynomial is equal to 

W?(s,-D m ) = C n (a,D m )P^ m) (s + er(D m )|D m ). (26) 

The partial polynomial part for the j-modified matrix D™ reads 

W?(s, Df ) = C n (a, Df )B«' m \s + cr(Df )|Df ), 

and similarly to the scalar case the partial vector Sylvester wave Wj n (s, D m ) for arbitrary j ^ 1 can 
be written as a linear superposition of the partial polynomial part for the matrix Dj" multiplied by 
the corresponding prime circulator: 

Pi-l 

Wj n (s,D m ) = W?(s - r • D^^Df )*j(s - r • D m ~^). (27) 
n=o 

Combining (|23|25() and (|27h we arrive at the final expression of the restricted vector partition 
function as a linear superposition of the vector Bernoulli polynomials of higher order multiplied by 
the vector prime circulators: 

T^(s,D m ) = ^ ^ Yl ^i n (s-r-D m -^J,D] n )*j(s-r-D m -^), (28) 

j n|=m— I ri=0 

where W r °(s,Dy i ) is given by 

W ^ D T) = nip^Dj") ^ ,m) ( 8 + <r WI D D- ( 29 ) 
It also can be written in the form of a mixture of partial waves 

W( S ,B m )= Y, C n (<x,-D m )W n (s,-D m ), (30) 

|n|=m— ( 
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where the superposition coefficients depend on the vector a and the partial wave is given by 

Pi-i 

W n (s,B m ) = J2Y1 7r _1 (p m ^ j )5n' m) (s + c(Df ) - r • D m ^|D] n )*j(s - r • D m ^). 
j n=o 

In the scalar case we have I = 1, D m = d m ,n = m — 1, and P m (a,dJ l ) = 7r(d™)a m . Thus, the 
expression (j29[) reduces to ()18[) . and noting that in the scalar case all pi = j, one finds that (|2T[) 
transforms into (|19p . 

We present several examples of application of the formula (|28p and show that it gives not only 
the vector partition function but also enables to find the chamber structure of the system. 



4.1 Example 1 

Consider computation of the vector partition function for the matrix (m = 3, 1 = 2) 



1 2 
1 1 



first using the partial fraction expansion suggested in [2]. An idea of the method is based on the 
definition of the generating function (I20j) . It is easily seen that VK(s,D) can be computed as the 
constant term in the expansion of the following expression: 



W(s, D m ) = const 



1 m 1 

-TT — 

1=1 



(31) 



Assuming all but one (say t{) components of the vector t to be constant expand the r.h.s. in (|3ip 
into the partial fractions in t\. This expansion contains both analytic and meromorphic parts w.r.t 
t\ = 0. The meromorphic part doesn't contribute to the ti-constant term, so it can be dropped. 
The constant term of the analytic part depends on the remaining components of t, so that such 
expansion eliminates t\. Applying this procedure repeatedly we eliminate all components of t and 
arrive to the result. 

Assume t\ constant and make partial fraction expansion w.r.t. obtaining 



1 



1 



1 



(l-t 2 )(l-iii 2 )i« 1-ti 
where MMP stands for the meromorphic part. Thus we have 

1 



VF(s,D- 



const/ 



const/ 



l - 1 2 i - ht 2 



l 



+ MMP, 



l - 1 2 i - ht 2 



Finding that for a > 
1 



const 



(1 - t 2 )(l - t)t a 
we obtain (note that si > 0) 



const 



■■> ^-?rConst t , 
(l-t?)(l-t!)tf 

1 - t S2+l 
(l-t 2 )(l-t)t^ 



l + 2a (-l) a 

+ .J - + - + mmp 



2(1 -t) 2 4(1 -t) 4(1 + t) 



a 3 + (-l) a 
2 + 4 ' 



const 



(1 - £ 2 )(1 - 



si 3 + (-l) gl 
2 4 



12 



and 

Jo, si - s 2 - 1 < 0, 

\ ^ + M T'" > -i-* a -i>o. 

Combining these results we arrive at the final expression for the partition function 



const* 



(1 — i 2 )(l - t)£si- S2 -i 



msD 3 )-/ t + ^^i Sl -s 2 -l<0, 

It is easy to check that the expressions in (|32p coincide for si = S2 and s\ = s 2 — 1. 

Turning to the formula (|30|) and having |n| = 1, we find that there are two partial waves 
corresponding to ni = (1,0), n 2 = (0,1). Computing the polynomial P3(q,D 3 ) = 2a\a 2 (ct\ + a 2 ) 
we obtain the coefficients: 

C (1 , 0) (a,D 3 ) = 01 C (0 , 1 )( Q ,D 3 )= a l 

Finding cr(D 3 ) = (3,2) and using (|26j) we have for the partial polynomial parts 

ttf' 0) (s,D 3 ) = C (1>0) («,D 3 ) L + , Wi°>%,l> 3 ) = C (0il) (a,D 3 )( S2 + 1). (33) 

In addition to the polynomial contributions there are nonzero terms corresponding to ji = (2, 1) 
and j 2 = (2, 2), while the term with j = (1, 2) doesn't contribute to the final result. 
The ji-modified matrix Dj^ reads (p3 = 2) 

D 3 - ( 2 ° 2 
" ^ o 1 2 

where the columns are sorted already. In order to use (|27[) we first compute the partial polynomials 
for the ji-modified matrix: 

wf' 0) (s,D 3 2jl) ) = q li0) («,D 3 ) (y + l) , <' 1} (s,D 3 2jl) ) = C (0il) (a,D 3 ) (f + ~) , 

and find 

i 

W ( n 2 ; i} (s,D 3 ) = Y,W?(( Sl -r,s 2 -r),nf 2A) )y 2 ( Sl -r), 

r=0 

arriving at 

wg^fcD 3 ) = ic (1 , 0) («,D 3 )^( Sl ), wg§(s,D 3 ) = iq 0il) (a,D 3 )* 2 ( Sl ). (34) 
The j2-modified matrix Dj^ reads (p 2 = P3 = 2) 

D 3 - ( 2 2 

U ( 2 ' 2 ) " V 2 2 

Repeating the computation we start with the partial polynomials for the j2-modified matrix: 



<' 0) (s,D 3 2j2) ) = C (li0) (a,D 3 ) (j + iy ^ 0,1) (s,D 3 2>2) ) = C7 (0 ,i)KD 3 ) (f + 
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The formula (|27p takes form 

i 

W^ 2) (s, D 3 ) = w?(( Sl - n, s 2 - n - r 2 ), D 3 2 2) )^ 2 ( Sl - n)* 2 (s2 -n- r 2 ), 

n ,r 2=0 

and we obtain 

w£$(s,D 3 ) = 0, wg$(s,D 3 ) = iq ,i)(« ) D 3 )* 2 ( Sl )M/ 2 ( S2 ). (35) 

Combining the expressions (133)) . (|34|) and (j35[) we arrive at the final result 

W(s,D 3 ) = C {li0) («,D 3 )W( 1 ' )(s,D 3 ) + q 0il) (a,D 3 )iy(°- 1 )(s,D 3 ), (36) 

where 

w^°>(s,d 3 ) = Sl + | + i* 2 ( Sl ), 

wW(s,D 3 ) = S2 + i + I^ 2 ( Sl ) + i^ 2 ( Sl - S2 ). 

Noting that v I / 2 (s) = (— l) s one can check by straightforward computation that the case a 2 /«i = 
produces the first line in (j32l) . while ai/a 2 = gives the last line. We see that proper choice of the 
vector a corresponds to the selection of the chamber. 

4.2 Example 2 

Consider another example for I = 2, m = 4 and 

D4= (l?Jl)' s = ( s ^). 
The straightforward computation using the partial fraction expansion gives the following result (see 

my- 



W{s,~D 4 



-1*2-^ + ^ + ^=-, *i/2-l<-a<-i + l, ( 37 ) 
^ + ¥ + 1, s 2 < si/2. 



v 2 1 2 

As |n| = 2, we have three partial waves corresponding to ni = (2,0),n 2 = (1,1), 113 = (0,2). 
Computing the polynomial Pi(a, D 4 ) = aia2(ai + a 2 )(2ai + a 2 ) we find the coefficients: 



q 2 ,o)(a,D 4 



2(ai + a 2 )(2«i + a 2 ) 



C ( i,i)KD)- (ai + a2)(2ai+a2r 

W^° 4 ) = 2(a l + a 2 a )(2a l + a 2 y 
Finding <x(D 4 ) = (4,3) and using ([26]) we have for the partial polynomial parts 

wf' 0) (s,D 4 ) = C (2iO) («,D 4 )^ 2 + 4 Sl + 0, 

wf' 1} (s,D 4 ) = C (lil) (a,D 4 )( SlS2 + ^i + 2 S2 + ^V 
ttf' 2) (s,D 4 ) = C (0 ,2)(«,D 4 )( S 2 + 3s2 + 2). 



(38) 
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In addition to the polynomial part only the term corresponding to j = (2, 1) produces nonzero 
contribution to the vector partition function, while two other terms with j = (1,2) and j = (2,2) 
don't contribute into the result. The j-modified matrix D™ reads (p^ = p\ = 2) 

D 4 _ ( 2 2 2 



( 2 .1) ~ V 1 1 2 

We first compute the partial polynomials for the j-modified matrix: 

■2 



wf'°>(s,D? 2il) ) = C (2 , 0) (a,D 4 )(^ + ^i + 2), 

^ J ( 8 . D W = ^(i,i) ( a » D ) ( — + Y + X + ¥ 

Wl°' 2 \ s ,-Df 2jl) ) = C (0 , 2) (c*,D 4 )(| + S2 + £), 



and use (|271) to have 



l l 



W ( n 2 ; i) (s,D 4 )= ^ ^^ ( ( Sl - r 3- r4 , S2 - r3 ),Dj 2il) )* 2 ( Sl -r3-r 4 ), 

r3=0 r-4=0 



and obtain 







;«,d 4 ; 


)*2(«1 




= 4^(1,1) 1 


;<*,d 4 ; 




<;5(s,d 4 ) 


= 0. 







(39) 



Combining the expressions (j38J) and (139p we arrive at the final result 

W(s,D 4 ) = C (2i0) (a,D 4 )^^^ 
with 

^( 2 <°)(s,D 4 ) = sl + i Sl + 7 - + ^ 2 ( Sl ), 

W^foD 4 ) = sis 2 + ^ + 2 S2 + ^ + ^ 2 (si), 
^(°' 2 )(s,D 4 ) = S 2 + 3s 2 + 2. 

We find that the case a%jct\ = produces the first line in (|37p . while «i/a 2 = corresponds to the 
last line. The second line is obtained as a real part of W(s, D 4 ) for a 2 /«i = — 1 ± i. Thus, we see 
again that the choice of the vector a corresponds to selection of the chamber. 



4.3 Example 3 

It should be noted that the formula (|28|) remains valid also in the case m = I. To illustrate it 
consider the case m = I = 2 with 




s = (si,s 2 ). 
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The solution of the corresponding linear problem reads x\ = S2, #2 = (si — S2)/2, thus the number 
of solutions equals to one for even nonnegative differences s± — S2 and zero otherwise. This can be 
written as 

V ^ 1 0, Si < 8 2 . K ' 

As m = I we have n = = (0,0), and all B^ n ' m \s\'DJ l ) = 1. The polynomial part is found as 

<(s,D 2 ) = °l ny (41) 

In addition to the polynomial part only the term corresponding to j = (2, 2) produces nonzero 
contribution to the vector partition function, while two other terms with j = (1,2) and j = (2, 1) 
don't contribute into the result. The j-modified matrix D™ reads (p = 2) 

^ , 2 2 ' 

and the corresponding term is found 



D f«> - V 2 



r=0 



Combining (|4ip and (|42p we arrive at the expression 

W(s,D 2 ) = °f + 
2(ai + a 2 ) 

This result also can be obtained from (|36|) applying to it the recursive relation (|21|) with C3 = (0, 1). 
It is easy to check that the choice 01/02 = corresponds to the first line in (|4*0|) . while 02/01 = 
produces the second one. 

4.4 Example 4 

Finally consider the case m = 4, Z = 3 with 

/ 2 1 \ 
D 4 = 112 , s = (si,s 2 ,s 3 ). 
\ 0010 / 

The straightforward computation using the partial fraction expansion gives the following result 

' (1 + (-l)-i+«+»s) (ttk^L + sx^ S2 > S3; Sl _ S2 + S3 _ 2 < 0, 

W(s,D 4 ) = (1 + ( _ 1)S1 + S2+S3) + , S2 > Sg , Sl _ S2 + S3 _ 2 > 0, (43) 

,0, s 2 < s 3 . 

There are three partial waves corresponding to ni = (1, 0, 0), 112 = (0, 1, 0), ni = (0, 0, 1). We drop 
the computation details and present the final result in the form 

^(s,D 4 ) = 31 , l . ■ j[3ai + 4a 2 + (-l) Sl («i + « 2 ) + (-1) S3 («2 + a 3 ) + 2a • s]. 

8(ai + a 2 )(«2 + a 3 ) 



The comparison of the real part of the above expression with (j43h leads to the following values of the 
vector a - the first line in (|4"3|) corresponds to 02/01 = 0,03/01 = ±i, the second one is obtained 
with 02/01 = ±2,03/01 = 1 =F h an d the last line is given by 03/01 = 03/02 = 0. 
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4.5 Discussion 

The formula (|28p provides an exact solution of the vector restricted partition problem through the 
vector Bernoulli polynomials of higher order. It also can be written in the form of a mixture of 
partial waves (|30f) and the superposition coefficients depend on the vector a. It appears that proper 
choice of this vector parameter produces the solution in one of the system chambers. The procedure 
for selection of the vector a will be published elsewhere. 

Knowledge of the vector a for each chamber gives way to determine chamber boundaries. The 
partition function (as well as its Sylvester waves) in two adjacent chambers should coincide at their 
boundary. Thus, setting the corresponding values to a and equating the real part of resulting 
expressions for W(s, D) one obtains the linear relations between components of the vector s which 
determines the location of the boundary. 

Thus, the formulas (|28h and (I30p appear to be the most general as they not only give the 
expression for the vector partition function in every chamber of the system, but also determine 
shape of each individual chamber. 

The derivation of (|28H was made in assumption that the matrix D is non-degenerate. The 
degeneracy of the matrix related to linear dependence of its columns doesn't affect the solution 
of the problem. It happens only if the matrix rows are linearly dependent. This dependence 
leads to effective reduction of the row number I and to additional linear conditions imposed on the 
components of the vector s. These conditions should be satisfied to have nonvanishing value of the 
vector partition function. This way the original problem is reduced to the non-degenerate case. 

References 

[1] M. Beck, Multidimensional Ehrhart reciprocity, 
J. Combin. Theory Ser. A 97 (2002), 187-194. 

[2] M. Beck, The Partial- Fractions Method for Counting Solutions to Integral Linear Systems, 
Journal of Discrete and Computational Geometry 32 (2004), 437-446. 

[3] M. Brion and Michele Vergne, Residue formulae, vector partition functions and lattice points 
in rational polytopes, 
J. Amer. Math. Soc. 10 (1997), 797-833. 

[4] L. Carlitz, Eulerian Numbers and Polynomials of Higher Order, 
Duke Mathematical Journal, 27 (1960), 401-423. 

[5] L. Carlitz, Multiplication formulas for generalized Bernoulli and Euler polynomials, 
Duke Mathematical Journal, 27 (1960), 537-545. 

[6] L. Dickson, History of the Theory of Numbers, V.2, Ch.3, 
Chelsea Publishing Co., NY, 1976. 
L. Comtet, Advanced Combinatorics, Ch.2, 
D. Reidel Publishing Co., Dordrecht, Holland, 1974. 

[7] F. G. Frobenius, Uber die Bernoullischen Zahlen und die Eulerischen Polynome, 

Sitzungsberichte der Koniglich PreuBischen Akademie der Wissenschaften zu Berlin (1910), 
809-847. 

[8] N. Nielsen, Traite Elementaire des Nombres de Bernoulli, 
Gauthier-Villars, Paris, 1923. 



17 



[9] N. E. Norlund, Memoire sur les Polynomes de Bernoulli, 
Acta Mathematica, 43 (1922), 121-196. 

[10] B.Y. Rubinstein, Expression for Restricted Partition Function through Bernoulli Polynomials, 
Ramanujan Journal, accepted. 

[11] B.Y. Rubinstein and L.F. Fel, Restricted Partition Function as Bernoulli and Eulerian Poly- 
nomials of Higher Order, 
Ramanujan Journal, 11 (2006), 331-346. 

[12] J. J. Sylvester, On Subinvariants, i.e. Semi-invariants to Binary Quantics of an Unlimited 
Order. With an Excursus on Rational Fractions and Partitions, 
American Journal of Mathematics 5 (1882), 79-136. 



18 



